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Abstract 



In this work, it is reviewed the cell model of adsorption by using the small system grand ensemble method. Under the 

^ c~| common assumption that the adsorbate phase is divided into identical and weakly-interacting subsystems, it is suggested 
Q-that the general multiparametric isotherm that arise from the theory may be used to fit the experimental data of adsorption 
Jh of gases and vapours on microporous adorbents (type I isotherm), even though these present heterogeneties like pore size 
Ch distribution. A simplified isotherm that reduces the number of adjustable parameters with respect to the general isotherm 
is proposed. Both isotherms, due to their relative high accuracy, can be used to estimate thermodynamical properties like 
(j isosteric and differential heats of adsorption. Also, a simple method is presented for systems that show an apparent variation 
^ in the coverage limit as function of temperature; for several systems, this method reduces the fitting deviations. Finally, several 
O applications to fitting data of adsorption, taken from literature, of some gases on activated carbon, molecular sieving carbon, 
c/3 silica gel, and pillared clays are presented. 
>> 
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1 Introduction 



C\| Correlating adsorption data obtained from experiments or 
,— 1 computer calculations is necessary to save time and efforts in 
additional experimentation or computing time. Because of the 
. complexity of equilibrium adsorption phenomenon, even to- 
1—1 day, the research of new models of adsorption is active. A 
good parameter-adjustable model must fit well the experi- 
mental data and correctly predict thermodynamic quantities; 
• • therefore, it is necessary to analyze the model and assess its 
. £^ capabilities. Many real- world applications of adsorption en- 
compass a range of saturation up to 10 6 and 10 12 in pressure 
H [1] . Therefore, in some cases, in order to fit the adsorption data, 



it may be necessary models that give experimentalist and pro- 
cess engineer the possibility to set the number of adjustable 
parameters. 

Many adsorbents used in industry and research present a 
certain degree of heterogeneity in terms of a pore size distri- 
bution and certain surface topography. Depending on each of 
these factors, different isotherm models are proposed. For 
example, lattice gas theories that simplify the structure of 
the adsorbent surface by taking into account the most im- 
portant adsorption sites of the adsorbent surface or adsorp- 
tion pores have been developed, and most of these theories 
are consistent with experiments and computer simulations 
[2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12]. A disadvantage of this kind 



of models is that many gas-solid systems present unique char- 
acteristics that must be considered in each particular model. 
Therefore, once a model has been developed for a certain gas- 
solid system it must be tested with adsorption and thermody- 
namic data. 

Another commonly used method to obtain isotherms for 
heterogeneous systems is the integration over a patchwise to- 
pography of adsorption sites [13]. Modified isotherms arise by 
assuming an adsorption energy distribution function and an 
ideal local isotherm like Langmuir. Nevertheless, although 
some of these equations work well for a large number of 
systems, they present limitations involved in using a par- 
ticular model, e.g., some models do not provide the correct 
Henry's law limit. Additionally, this method is somewhat 
difficult because the complexity of the adsorption energy dis- 
tribution makes it difficult to obtain an analytical expression 
for the transformed isotherm. The Sips[14], Toth[15], and 
UNILAN[16] isotherms are the most widely used isotherms 
of this type. 

Recently, a method proposed by Reiss and Merry [17] to 
study the small system grand ensemble (SSGE)[17, 18, 19, 20] 
has been useful to analyze the properties of these small subsys- 
tems that all together form a homogeneous phase. Initially, this 
model was proposed to study rigid spheres, but, in principle, 
it can be applied to other systems [19], as shown in this work 
for the adsorption case. The advantage of this SSGE method is 
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that it can be studied the implications of assuming interactions 
between a particular subsystem and its surroundings, allow- 
ing us to analyze in detail the properties of these small systems 
and their consequences for the adsorption of gases. 

In this work, it is briefly reviewed the statistical mechanical 
cell model of adsorption[21, 22, 23] by using the SSGE method. 
The general isotherm that arises from the cell model has been 
useful to fit experimental data of adsorption of gases on ze- 
olites [24, 25, 26, 27, 28, 29, 30, 31]. Here, an accurate and 
simplified isotherm is proposed, this isotherm is similar to the 
empirical isotherm proposed by Ruthven [32]. It is shown that 
both the simplified and general isotherms might be applied to 
the adsorption of gases on heterogeneous microporous adsor- 
bents. To prove this hypothesis, several experimental adsorp- 
tion isotherms are correctly fitted by using the isotherm equa- 
tions studied here. From the fitting results, it is suggested that 
the isotherm equations treated here can be used to fit type I ex- 
perimental adsorption isotherms and predict thermodynamic 
properties. Furthermore, extrapolations to temperatures near 
the experimental temperature range are possible. Because of 
the complexity of heterogeneous adsorbents, it is argued that 
the method presented here is semi-empirical and the cells en- 
semble can be considered as an abstract and idealized grand 
ensemble of small subsystems. 



2 Theory 

2.1 The General Adsorption Isotherm 

Consider a single-component gas in equilibrium with an adsor- 
bate phase. Suppose that the adsorbate phase can be divided 
into M identical cells that weakly interact between each other. 
The total number of cells is temperature-independent and the 
cell volume as well. Now, let us regard a particular cell as a 
subsystem of interest and the remaining M - 1 cells that sur- 
round this subsystem as a bath. It is assumed that the effect of 
the adsorbent on the adsorbate manifests only by means of an 
interaction potential and the adsorbate presence on the adsor- 
bent surface (or cavities) does not affect its properties. Define 
the following partition functions: 
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where A is the deBroglie thermal wavelength, V and N are 
the adsorbate phase volume (or area) and the total number of 
adsorbed molecules, respectively, v s and n are the volume (or 
area) and number of molecules of the subsystem, respectively, 
and d is the adsorbate phase dimension. The potential energy 



may be decomposed into three contributions: 

^(R N ) = 36(Ri, . . . , R„) + <%(R„ + i, . . . , R N ) 



(3) 



the term is the potential energy of a subsystem and % 
is the bath potential energy, these potentials are the sum of 
the molecule-molecule and the molecule-solid interaction po- 
tentials. The term represents the interaction between the 
subsystem and the bath. Consider the usual probability den- 
sity for the canonical ensemble [23, 33]: 



p(P N ,R N ) = Cxexp[- 



(4) 



where C is a constant. If we integrate over the subsystem's 
phase space and the bath's phase space, we obtain the proba- 
bility of finding n molecules in the subsystem volume v s and 
N - n molecules in the bath volume V - v s : 



p n =CxQ*(T,V-v s ,N-n,v s ,n) 



(5) 



It is plausible to assume that a maximum of n s molecules can 
adsorb on the subsystem (the adsorbent has a saturation limit). 
Therefore, in order to eliminate the arbitrary constant C, we 
can write: 

p. = ») 

Pj may be written as follows: 

Q*(T,V-v S/ N-j,v Sr j) 



Pi =C x 



Q{T,V-v s ,N-j)Q{T,v s ,j) 
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By invoking the chemical potential definition we obtain: 



Q(T, V -v s ,N - j) = Q(T, V - v s , N) exp (i^) 



By defining 
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it is obtained 

Pj = CxQ(T,V-v s ,N)exp 
where A = e^'. Hence, we have: 
exp 



A' 



(7) 

(8) 

(9) 
(10) 
(11) 
(12) 

(13) 



n^ s -A° n -A s n 



knT 



1 + Yj eX P 

7=1 



■A a -A s 

1 1 



(14) 



2 



where ju° g is the gas-phase chemical potential at the reference 
pressure P°, and a is the activity: 



a = 



J_ 



The term A" may be regarded as a free energy of interaction 
between the subsystem and the bath [17]. If we assume that 
the bath and the subsystem weakly interact, then is small. 
Now, we define: 



Kj{T) = exp[ 



(15) 



If the system may be divided into identical weakly-interacting 
[34] subsystems then we can state that H ~ <£ M where [22]: 



£ = 1 + Ki,l 
7=1 



(16) 



here, £ is the subsystem's grand canonical partition function 
and S is the grand canonical partition function. We can calcu- 
late the saturation (q) as follows: 



7=1 
1 + 

7=1 



(17) 



The term Kj can be considered as an adsorption equilibrium 
constant. Therefore, by similarity with chemical equilibrium 
constants, we can propose the following relation: 



h\K° 



Ahj 
RT 



(18) 



where h\K° and Ahj are a change of entropy and enthalpy, 
respectively, related with the adsorption of ; molecules on a 
representative microscopic subsystem at the reference pressure 
P°. Eq. (17) was proposed by Langmuir [35] for the case in 
which a site can hold several molecules and the adsorbent is a 
set of non-interacting sites. In light of his proposal, the term 
Kj is partitioned as follows: 



K 



1 i 

1 1=1 



In a similar fashion, Rj behaves as Kj: 

R ; = R°exp(— I) 



where 



jK° 

R° = — - 

I K° 
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Ahj = Ahj - Ahj-i 



(19) 

(20) 

(21) 
(22) 



Thus, in order to fit an individual isotherm using Eq. (17), the 
adjustable parameters are {Rj} instead of {Kj}. For the case in 
which several isotherms at different temperatures have been 
measured, the adjustable parameters are {lnR°}, {Ahj}, and q mr 
a total of 2n s + 1 parameters. 

The present model suggests that £ — > H as the subsystem 
volume is increased. Thus, in order to predict q, we should 
assume that the maximum number of molecules n s is large 
enough such that A a . — > for all j [17]. However, if we as- 
sume that each Rj is an adjustable parameter, this would give 
us a high number of adjustable parameters with poor statis- 
tical confidence [31]. A solution for this problem could be 
assume n s as a small number of molecules that adsorb on 
a microscopic imaginary subsystem and the parameters {Rj} 
as representative parameters of the experimental adsorption 
isotherm. Thus, we are assuming that the adsorbate system 
can be divided into an abstract and idealized grand ensemble 
of subsystems. The characteristics of the adsorbate+adsorbent 
system like molecule-surface and molecule-molecule interac- 
tions are included abstractly in these subsystems. This is an 
ad hoc guess because we are not claiming that this idealized 
system exists, but this assumption is a convenient picture that 
allows us to apply the cell model even to heterogeneous sys- 
tems despite this model was originally conceived to study 
relatively homogeneous systems like zeolites. 

Accordingly with Eq. (6), the probability that a subsystem 
with n adsorbed molecules will be found is (assuming ideal 
gas phase): 

„ _ k„(p/pT (23) 



Eq. (17) can be written as: 



S 7=1 7=1 



(24) 



here, 6 = q/q m , and /; is a fraction of molecules that can be 
found in subsystems with j molecules. It is known that, at 
low pressures, the leading term in Eq. (24) is f\ (as stated by 
the Henry's law), and, as the pressure increases, the leading 
term is fi. Similarly, in a certain pressure range, each /; term 
significantly contributes to the fractional coverage. Therefore, 
each parameter Kj is estimated with data taken in a certain 
range of the saturation q. Hence, for intervals where there is not 
enough data, these parameters will give large uncertainties. 

The general isotherm (Eq. (17)) and its simplifications have 
been useful to fit sets of isotherms of adsorption of gases 
and vapors on zeolite adsorbents at different temperatures, 
where the experimental isotherms show clearly a saturation 
limit, or are reported within the same range of saturation 
[24, 25, 26, 27, 28, 29, 30]. However, a great number of adsorp- 
tion systems do not show this behavior because experiments 
typically are performed within the same pressure range; there- 
fore, the isotherms apparently show that the saturation limit 
depends on temperature (i.e., q m tend to increase with a de- 
crease in temperature). Now, consider a certain isotherm at 
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temperature 7\, if we consider that n s = 3, we have the follow- 
ing adjustable parameters: K\, K2, K3, andq m (Ti). On the other 
hand, suppose that we have a second isotherm at temperature 
T2 and its apparent saturation limit is q m (T2)', if Ti < T\, then 
qm(T2) > cjm(Ti). This would suggest that, as the temperature 
decreases, new subsystems are created, but this would be a 
physical inconsistency. To solve this problem, we introduce a 
new equilibrium parameter K, h+ \ that take into account the ad- 
sorption of an extra molecule on a subsystem at high pressures 
and low temperatures. Because this new parameter describes 
the adsorption at these conditions, it is possible that there is 
not enough experimental data to estimate lnR° , and Ah„ s+ \, 
then it is plausible to assume that R, k+ i = R, k and write 



K 



n s +l — 



n s + 1 



(25) 



Thus, we are assuming that the change of both entropy and 
enthalpy of adsorption of the (n s + l)th molecule on a micro- 
scopic cell is the same for the adsorption of the n s th molecule 
(i.e., RlnR° = RlnR° s and A/z„ s+1 = A/x„ s ). This result is 
generalized as follows: 



Now, E becomes: 



and the isotherm is: 
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2.2 A Simplified Adsorption Isotherm 

In order to reduce the number of adjustable parameters and 
reduce confidence intervals around parameters estimates, let 
us re-express Eq. (17) as follows (assuming ideal gas phase): 
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here a = P/P°. It can be shown that: 



fc B T 2 --ln 



Z{T,v s ,j) 



dT [(k B T/P°)i 



where 



= <^s(/)> - /for 



(29) 



foT 



(30) 



because the adsorption energy (^s(j)) is usually greater than 
/foT, then we can assume that the latter term is a constant /foT, 
where T is a mean temperature. Integration of Eq. (29) gives: 



In 



Z(T,v s ,j) 



(kuT/P )' 



= K° exp 



(31) 



As an special case, let us assume that the potential energy (;') 
can be written as follows: 



= Yj ^f-m(R«) + Yj ^m-m(R„ Rm) (32) 

n=l n>m 

where ^f-m and ^f-m are the field-molecule and the molecule- 
molecule interaction potentials, respectively. Now, we obtain 

ma)) = im-u) + ^-t^Vm-m) (33) 

For the sake of simplicity let us assume that (^-m) and (^m_m) 
do not depend on j and there is a temperature T° such that 
K° = (K°y, where K° is the pre-exponential term in the Henry's 
constant. Under these assumptions, it is obtained the adsorp- 
tion isotherm: 



£- -7-r-; 
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1 + £ (K°a)> 
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(34) 



where 



E M _ F = -N A m- M ) 
Em-m = -A^a('^m-m) 



(35) 
(36) 



here, Na is the Avogadro's constant. If (%i_m) and (^m-m) 
depend on /, we can introduce the equation: 



(K°a)/ 



M s + / -Ahj/RT 



(37) 
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comparing with Eq. (27) we have: 



K j = 



For this simplified isotherm, from Eq. (26), we have: 

^o e -(A/ I „ s -AV_ 1 )/RT 



« s +/ — 



n s + / 



-K, 



n s +/-l 



(38) 



(39) 



Here, the fitting parameters are {Ahj), K°, and q m , a total of 
n s +2 adjustable parameters. We can connect this result with 
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the empirical isotherm proposed by Ruthven [32] (1 = 0): 



It is easy to show that [37, 38] 



1 = 



qmM_____ 

l+ Li {KP) 



Ai =1 



(40) 



in this Eq., K is the Henry's constant, and Aj is an empirical 
parameter. The parameter Aj is expressed as follows: 



, Nhj, 



and 



K = 



K° 
P0 



(41) 



(42) 



thus, the Henry's constant in Ruthven's Eq. is replaced by the 
pre-exponential term K°/P°. 

If the term k^T is not precise, Eq. (37) can be written as 

*+» e -Au,/RT R o ■ 

2-i (/' - 1)! v T' ' 

g=-*4 ^ (43) 

«s + ' >%+! a -Auj/RT go a ■ 



T 



where 

r - p ( 44 ) 

Any = N A (®s(f)) (45) 
here, T° is a reference temperature (e.g., 1 K). 

2.3 Isosteric Heat of Adsorption 

From Eq. (15) we have 



(46) 



where Uj and hi are the energy of j molecules in a subsystem 
and the molar enthalpy of the gas phase at P°, respectively, 
Ahj is the change of enthalpy related to the adsorption of j 
molecules on a subsystem. It can be noticed that, in princi- 
ple, Ahj is temperature-dependent, but for the sake of sim- 
plicity, we have assumed that both h® and Uj are temperature- 
independent. The isosteric heat of adsorption can be calculated 
by means of the following formula [36]: 



q st = h g -u a 



(47) 



where he and w a are the molar enthalpy of the gas phase and 
adsorbate molar internal energy respectively. For any function 
gj, define: 



(jAhj) - n s e(Ahj) 
{f)o-n 2 s d 2 



(48) 



where 6 = q/q m . From the above eq, the two following limits 
arise: 



limfl st 

9^0 



lima 



-A/ti 

-Ah ns 



(49) 
(50) 



Equivalently, if Eq. (43) is used, the differential heat of adsorp- 
tion can be calculated as follows: 



(jAuj)o - n s 6(Auj) 
{j 2 )o-n 2 s d 2 



(51) 



For the experimental isotherms studied here, both simplified 
isotherms (Eqs. (37) and (43)) give the same results. Here- 
inafter, Eq. (37) will be used due to its simple relationship 
withEq. (48). 

3 Results and Discussion 

To fit the experimental data cited here, the unweighted least- 
squares Levenberg-Marquardt algorithm was used. Details of 
the results of fitting Eqs. (37) and (27) to several experimental 
adsorption data are shown in Tables 1 and 2; the parameters 
are tabulated with their corresponding marginal confidence in- 
tervals [39, 40] . The study of these confidence intervals is nec- 
essary to avoid over-interpretation of the parameters [41] and 
to estimate uncertainties in predicted thermodynamic proper- 
ties. The parameter li\K° (or lnX°) was used instead of K° (or 
K°) due to the restriction K° > (for all j). To analyze the fit- 
tings precision, the following deviation percentage parameter 
was used: 



D 



100% 
N T 



Ntemp N P (j) 

EE 
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q{Pi,Tj)-q exp (P ir Tj) 



qexp(Pi/ Tj) 



(52) 



where Nt em p is the number of isotherms, Np(j) is the number of 
experimental data taken at Tj temperature, Nj is the total num- 
ber of experimental data, q(Pi, Tj) is the calculated saturation, 
and q ex p(Pi, Tj) is the experimentally measured saturation. It 
can be noticed in Table 1 that the deviation parameter D in all 
cases is less than 5%. 

In Fig. 1 the fitted and experimental isotherms of SF(, ad- 
sorption on pillared clay (designated as W-A(673) [42]) are 
shown 1 , it was used Eq. (37). It can be noticed a good agree- 
ment between the experimental data and Eq. (37). For this set 
of isotherms, the correction term (Eq. (39)) was used (1 = 1) 
and it was obtained a deviation of 2.97 %. If this correction 
is not applied (/ = 0) the deviation is 5.2 %. Therefore, in 



^or this system, the saturation is expressed as a volume at 101.325 kPa 
and 273.15 K 
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Table 1: Estimated parameters for several systems 

q m (orv) lnR° -10 -3 x Ahj/R 


using 


Eq. 37. aM 




D 


Temperature 


Pressure 


System 


mmol g 1 (or cm 3 g 1 ) 


K 


% 


range (K) 


range (kPa) 


SF 6 +W-A(673)[42] 


23.89 ± 0.44 -13.67 ± 0.16 


3.107 ± 0.055 


2.97 


266.5-297.5 


0.1-100 


1 = 1 




6.112 ± 0.077 
8.08 ± 0.50 
11.62 ± 0.32 








C 3 H 8 +W-A[42] 


34.6 ± 3.9 -14.65 ± 0.30 


3.391 ± 0.096 


3.2 


267-298 


0.1-100 


1 = 1 




7.472 ± 0.097 
10.8 ± 0.23 

14.8 ± 0.60 

17.9 ± 0.48 








HFC-227ea+AC[43] 


3.87 ± 0.27 -16.67 ± 0.27 


6.066 ± 0.098 


2.36 


283.15-363.15 


0.01-100 


1 = 




11.39 ±0.19 
16.38 ± 0.32 
20.82 ± 0.45 
24.71 ± 0.68 








HFP+AC[43] 


3.87 ± 0.20 -15.66 ± 0.20 


5.429 ± 0.075 


2.53 


283.15-363.15 


0.01-100 


1 = 




10.15 ±0.15 

14.68 ± 0.23 

18.69 ± 0.32 
22.33 ± 0.45 








CO z +ZC[44] 


4.82 ± 0.39 -15.71 ± 0.29 


5.03 ± 0.10 


3.84 


273.15-353.15 


0.05-100 


1 = 1 




9.39 ± 0.20 
13.45 ± 0.31 
17.00 ± 0.43 








C0 2 +SG[45] 


14.5 ±1.0 -15.22 ±0.10 


3.068 ± 0.030 


1.47 


278-328 


50-3400 


1 = 4 




5.759 ± 0.052 
8.454 ± 0.094 
10.87 ±0.12 









" In all cases, P° is expressed in kPa. 

" Abbreviations: activated carbon (AC), 1,1,1,2,3,3,3-heptafluoropropane (HFC-227a), hexafluoropropene (HFP), zeocarbon (ZC), silica gel (SG). 
c The parameters Ah, are tabulated in increasing order of /', e.g., for the SF(j+W-A(673) system, A/ij/R = -3.107 x 10 3 , A/i2/R = -6.112 x 10 3 , and so on. 
d For SF6+W-A(673) and C3H8+W-A systems, the saturation is expressed in cm 3 g _1 at 101.325 kPa and 273.15 K, and for the remaining systems it is 
expressed in mmol g _1 . 



Table 2: Estimated parameters for the adsorption data presented by Watson et al. [46] using Eq. 27 (1=0). " ^ 







lnR° 

; 


-10" 3 x Ahj/R 


D 


Temperature 


Pressure 


System 


mmol ■ g _1 




K 


% 


range (K) 


range (kPa) 


CH 4 +MSC 


4.412 ± 0.067 


-14.4 ± 1.5 
-13.88 ± 0.48 


2.81 ± 0.41 
1.92 ± 0.13 


2.29 


148-298 


1-4000 


CO2+MSC 


5.87 ±0.14 


-13.7 ±1.3 
-14.00 ± 0.83 
-16.07 ± 0.82 
-19.46 ± 0.58 


3.26 ± 0.40 
2.94 ± 0.25 
3.29 ± 0.22 
3.58 ± 0.17 


1.41 


223-323 


25-5200 


N2+MSC 


5.75 ± 0.12 


-15.6 ±1.1 
-14.03 ± 0.45 
-18.3 ± 1.0 


2.71 ± 0.31 
1.725 ± 0.099 
1.79 ±0.15 


3.14 


115-298 


0.01-5000 



" In all cases, P° is expressed in kPa. 

" Abbreviation: molecular sieving carbon (MSC). 

c As in 1, the parameters lnR° and Ahj are tabulated in increasing order of ;', for example, for the CH4+MSC system, lnR° = -14.4, lnR 2 = -13.88, 
Ahi/R = -2.81 x 10 3 , Ah 2 /R = -1.92 x 10 3 . 
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Figure 1: Comparison between the experimental data of SF6 
adsorption on W-A(673) [42] and Eq. (17); symbols: experi- 
ment, solid line: Eq. (37). 



this case, the correction term gives a better result. Fig. 1 sug- 
gests that Eq. (37) is useful for isotherms with complex shapes 
like those reported in Ref. [42]. The experimental isotherms 
shown in Fig. 1 were also fitted using the 5-parameter Toth 
isotherm and it was obtained a deviation of 17.8%; a value 
that is grater than that obtained with Eq. (37). On the other 
hand, Eq. (27) could be also used to fit these experimental 
isotherms, but more parameters would be necessary and this 
would lead to large error bars and thermodynamic properties 
with large uncertainties. Bandosz et al. [42] reported ad- 
ditional single-component experimental isotherms of propane 
and sulfur hexafluoride adsorption on various heat-treated pil- 
lared clays; a total of six adsorbate+adsorbent systems were 
studied. Given that the results are similar, here only two ana- 
lyzed systems are reported in Table 1 . But, for each of the six 
adsorbate+adsorbent systems the deviation was within 1 and 
3%. 

Fig. 2 shows the estimated isosteric heat of adsorption as a 
function of adsorbed volume and the corresponding marginal 
confidence intervals [39, 40] for the SF6+W-A(673) system. This 
Fig. shows that the isosteric heat of adsorption is nearly con- 
stant around 25.5 kj/mol. Bandosz et al. [42] used a virial type 
isotherm [47, 48] and the adsorption isosteric heat vs adsorbed 
volume plot obtained by them presents some oscillations near 
zero saturation around 25 kj/mol; this result does not agree 
with that obtained here by means of Eq. (48). These differences 
could be due to the method used by Bandosz et al. [42], they 
fitted small subsets of 15 adsorption data by using the virial 
type isotherm and used four parameters to fit each subset. Al- 
though this method is appropriate, the error bars reported by 
them correspond to standard deviations in calculated isosteric 
heats of adsorption. If 95 % marginal confidence intervals are 
used, the error bars at low coverage would be even larger than 
those reported by them and the isosteric heats reported here 
would fall within such uncertainty region. 
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Figure 2: Plot of isosteric heat of adsorption vs adsorbed vol- 
ume, the system temperature is 283.0 K; error bars correspond 
to 95 % confidence intervals. 
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Figure 3: Isosteric heat of adsorption for SF6+WA-(673) as a 
function of temperature at three different saturations. 



To illustrate the dependence of the isosteric heat of adsorp- 
tion in terms of temperature, plots of ^ st vs T at three different 
saturations are shown in Fig. 3. As expected, the isosteric heat 
of adsorption does not significantly change within the exper- 
imental range of temperatures. This is a typical behavior of 
the isosteric heat of adsorption that is observed in both exper- 
iments and molecular simulations. In the case in which it is 
necessary to take into account the dependency of isosteric heat 
of adsorption on temperature, the well known thermochem- 
ical formula [49] to calculate chemical equilibrium constants 
could be used to estimate In Kj as a function of temperature in 
Eq. (15); an empirical model for the heat capacity would be 
necessary. 

The fitting of the experimental isotherms of adsorption of 
1,1,1,2,3,3,3-heptafluoropropane (HFC-227ea) on activated car- 
bon is shown in Fig. 4. As in the previous case, there is a good 
agreement between model and experiment; the correction term 
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was not necessary (/ = 0). Due to the complex pore structure 
of activated carbon, it is difficult to consider a subsystem as 
any specific region of the real adsorbent, for this reason it is 
convenient to regard each cell as a representative subsystem. 
The isosteric heat of adsorption calculated by using Eq. (48) 
and Toth equation are shown in Fig. 5. It can be noticed that 
both curves fairly agree. The differences are attributed to the 
fact that Yun et al. [43] only used two isotherms to calculate 
the isosteric heat of adsorption; they employed the parame- 
ters of each isotherm and the Clausius-Clapeyron equation to 
obtain the q st vs q plot, in contrast, here the complete set of ex- 
perimental isotherms was used. A problem with the method 
used by Yun et al. [43] is that the isosteric heat of adsorption 
is influenced by uncertainties in adsorption data, and hence, 
very reliable isotherms are necessary to determine q st . 

In Fig. 6, it is shown the adsorption isotherms of CO2 on zeo- 
carbon [44] and the fractions of molecules distributed among 
cells at 273.15 K and 313.15 K. The zeocarbon synthesized by 
Lee et al. [44] is a zeolite X/activated carbon mixture composed 
by 38.5 mass % zeolite X, 35 mass % activated carbon, 10 mass 
% inert silica, and 16.5 mass % zeolite A and P. For this system, 
Eq. (37) accurately fits the experimental data (D = 3.8%), and 
the correction term (/ = 1) in this case reduces the deviation 
by 2 % with respect to / = 0. The deviation is less than that 
obtained by using the Toth equation, which is one of the most 
used equations to fit these kind of adsorption data. Given that 
the zeocarbon is a zeolite/activated carbon composite, from 
these results, it can be inferred that the method used here is 
empirical. Consequently, the assumption that the adsorbate 
phase may be divided into identical weakly-interacting sub- 
systems and the fictitious representative grand ensemble of 
subsystems are just convenient pictures for fitting and correla- 
tion purposes. 

As it was mentioned, at low pressures the leading term is fa 
and each fa significantly contributes in a certain region of the 
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Figure 4: Comparison between the experimental data of HFC- 
227ea adsorption on activated carbon [43] and Eq. (37); sym- 
bols: experiment, solid line: Eq. (37). 
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Figure 5: Isosteric heat of HFC-227ea adsorption on activated 
carbon calculated; dotted/dashed line: calculated by Yun et al. 
[43], solid line: Eq. (48), dotted lines: 95 % marginal confidence 
intervals for Eq. (48). 



isotherm. It can be noticed in Fig. 6a that near 100 kPa, the satu- 
ration (q) changes in approximately 0.5 mol kg -1 . Accordingly 
with the fitting method used here, at the lowest experimental 
temperature, for which at high pressures the saturation is ap- 
proximately q m , the leading term is / Hs +i, as shown in Fig. 6b. 
However, as temperature increases in the high experimental 
pressure region, the isotherm is a combination of several frac- 
tions fa. For example, Fig. 6b shows that the isotherm at 313.15 
K in the high pressure region is a combination of fa, fa and fa. 
When temperature is further increased, the term fa does not 
contribute to the isotherm. If the correction term is not used 
and n s = 4, the plot of fa vs P is quite similar to that shown 
in Fig. 6c for fa, and thus, this term (fa) does not significantly 
contribute to the isotherm at 313.15 K. For this reason, it was 
necessary to propose the correction terms shown in Eqs. (26) 
and (39), these Eqs. assure that high-order fractions depend 
on parameters that substantially contribute to estimate each 
isotherm within the experimental temperature range. Because 
of the specific characteristics of each adsorbate+adsorbent sys- 
tem, it is difficult to establish a priori whether the correction 
term is necessary, it must be tested whether this correction 
reduces the fitting standard deviation. 

In Table 2, the results of the fittings of data reported in Ref . 
[46] are presented. In this case, Eq. (27) was used without 
any correction. For the CO2+MSC and N2+MSC systems, Eq. 
(27) gives better results (in terms of deviation) than the sim- 
plified model (Eq. (37)). With Eq. (37) the deviation is around 
3.5% (/ = 0), whereas Eq. (17) gives deviations less than 3 %. 
Also, the deviation obtained using Eq. (27) is less than that 
obtained by using the Toth equation; this is also confirmed by 
F- tests. The confidence intervals for lnR° are larger than those 
obtained for lnX° in Table 1, this could be caused by parame- 
ter correlation effects. To fit this set of isotherms, Watson et al. 
[46] used the Toth isotherm and also obtained large confidence 
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Figure 6: a) Comparison between the experimental data of 
CO2 adsorption on zeocarbon [44] and Eq. (17); b) distribution 
of molecules among cells at 273.15 K; c) same as b) for 313.15 
K. 



intervals for the K° parameter. For the CH4 system, the Toth 
isotherm gives better results than Eqs. (17) and (37), this is due 
to an isotherm at 148 K that present an apparent saturation 
limit that is less than each of the apparent saturation limits 
of the isotherms at temperatures greater than 148 K. If this 
isotherm at 148 K is eliminated, the results obtained by using 
the 5-parameter Toth Eq. and Eq. (27) are quite similar. 



The advantage of the present models is the flexibility of 
setting the number of adjustable parameters; this condition is 
essential to fit isotherms with complex shapes. It is known that 
some of the most widely used empirical isotherms to describe 
type I isotherms are Sips, Toth, and Dubinin type [50, 51, 52] 
isotherms. These Eqs. do not reduce to the correct Henry's 
law limit; except the Toth isotherm, but this overestimates the 
Henry's constant [53] . In contrast, the isotherm models studied 
here have the advantage that they present the correct Henry's 
law limit and thus they can be used to estimate this constant 2 . 
Despite its advantages, the model studied here cannot give 
site energy and pore size distribution. However, although 
the model does not explicitly take into account the adsorbent 
heterogeneity, the number of parameters and the degree of 
the polynomial £, which are related to the subsystem size, 
give an idea of the variety of adsorbent sites and molecule- 
molecule interactions because the system size is related to these 
characteristics of the adsorbate+adsorbent system. 



4 Conclusions 

The cell model of adsorption was reviewed and analyzed by 
using the SSGE method and a simplified accurate isotherm 
model is proposed (Eq. (37)). It was found that the Eqs. (27) 
and (37) can be applied to fit experimental data of adsorption of 
gases and vapors on microporous heterogeneous adsorbents; 
these isotherms give relatively accurate fittings. A simple cor- 
rection that improve the fitting results is proposed. However, 
in some cases this correction may not be necessary, is must be 
tested. Additionally, for systems in which the experimental 
temperature range is large, it is suggested that the dependence 
of Ahj on temperature must be considered and a model for 
both gas and adsorbate phase heat capacity could be required. 

Finally, the model based on an idealized and abstract grand 
ensemble of subsystems presented here is useful to correlate 
single-component adsorption data and the corresponding ther- 
modynamic properties of interest. Specifically, the model is 
useful to fit type I isotherms, i.e., the reversible micropore 
filling region of the experimental adsorption isotherm. The 
advantages of the present model are the correlation accuracy 
of both adsorption and thermodynamic data, the flexibility to 
set the number of adjustable parameters and to consider varia- 
tions of In Kj with temperature, and the possibility of regarding 
adsorptive as a real gas phase. However, the method studied 
here does not consider explicitly the pore size and adsorption 
energy distribution, but the size of a representative subsystem 
gives an idea of the adsorbent heterogeneity because the size 
of the subsystem depends on this factor. 
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